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SUMMARY 

The trajectories of water droplets in a compressible-air flow field 
around a cylinder were computed with a mechanical analog. The results of 
the calculations at approximately the flight critical Mach number were 
compared with calculations of trajectories in an incompressible flow 
field. For a cylinder, the effect of compressibility of the air on the 
droplet trajectories was negligible up to the flight critical Mach number. 
The results obtained with the cylinder were extended to airfoils. This 
extension is possible because the incompressible flow fields of both 
cylinders and airfoils are similarly altered by compressibility. 


INTRODUCTION 

The analytical study of ice formation on airfoils and other aero- 
dynamic bodies exposed to icing clouds requires the calculation of the 
trajectories of the droplets with respect to the body. Previous inves- 
tigators have calculated the water -droplet trajectories with respect to 
cylinders (references 1 to 3) and with respect to airfoils (refer- 
ences 4 to- 6). All the calculations in the references cited were made 
with the assumption that the air surrounding the aerodynamic bodies 
behaved as an incompressible fluid. An evaluation of the effect of 
the compressibility of air- in the flow field on the trajectories of 
cloud droplets is important, because many airfoils are now operated at 
high speed. 

A knowledge of the air flow around the aerodynamic body is required 
in the calculations of the droplet trajectories, because the motion of 
a droplet is governed by the drag forces imposed on the droplet by its 
motion relative to the air moving along the air streamlines. The droplet 
momentum tends to keep the droplet moving in a straight path, while the 
drag forces tend to force the droplet to follow the streamlines. 
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Two-dimensional incompressible flow fields may be determined either 
from known mathematical relations or from electrolytic tank methods for 
most aerodynamic bodies. The compressible-air flow fields are extremely 
difficult to obtain for most aerodynamic bodies and are usually approxi- 
mations or incompressible flow fields with corrections applied. A 
reasonably simple method for obtaining droplet- trajectory results in a 
compressible -air flow field is to evaluate the compressibility correc- 
tion to be applied to trajectory data calculated for an incompressible 
flow field. Although the compressible-air flow field around a cylinder, 
which is obtained by applying correction terms to the incompressible 
flow field, is laborious to calculate, the cylinder flow field provides 
a relatively simple method for evaluating the correction required for 
compressibility. 

In order to evaluate the magnitude of the compressibility effect 
on droplet impingement, trajectories in a compressible flow field around 
a cylinder were calculated at the NACA Lewis laboratory and the results 
compared with those of trajectories calculated in an incompressible flow 
field (reference 3). The compressible -air flow field was determined 
only for the cylinder flight critical Mach number , which is defined as 
that airspeed which results in sonic velocity at some location on the 
cylinder . 

The results of the study of the effect of the compressibility of 
air around a cylinder are extended to include airfoils. The extension 
is based on a comparison of the compressible -air flow field around a 
cylinder at the cylinder flight critical Mach number with the 
compressible -air flow field around an airfoil at the airfoil critical 
Mach number. 


SIMBOLS 

The following symbols are used in this report: 
a droplet radius, ft (3,048X10^ microns) 

C-q drag coefficient for droplets, dimensionless 

E collection efficiency based on cylinder radius, dimensionless 

2 0 ^% 

K inertia parameter, ^ ^ , dimensionless 

L cylinder radius or airfoil chord length, ft 
M Mach number, dimensionless 
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Re 


Re 0 


r 


t 

U 

u 


v 


V 




a 

r 

9 

P 

P 

T 


cp 


local Reynolds number with respect to droplet, 2a p a v/p, dimension- 
less 


free -stream Reynolds number with respect to droplet, 2aU 
dimensionless 



radial distance from cylinder axis to point in flow field, ratio to 
cylinder radius 

airfoil thickness, ratio to chord length 

free -stream velocity, ft/ sec 

local air velocity, ratio of the actual velocity to the free-stream 
velocity 

local droplet velocity, ratio of actual velocity to free-stream 
velocity 

local vector difference between velocity of droplet and velocity 
of fluid, ft/ sec 

rectangular coordinates, ratio of actual distance to cylinder 
radius L or airfoil chord length 

angle between x-axis and radial line to point in flow field, deg 


ratio of specific heats 


impingement angle on cylinder, deg 


viscosity of air, slugs/(ft) (sec) 
density, slugs/cu ft 

dimensionless time scale, ratio to time required for droplet to 
travel distance L at velocity U 


Re 2 
ne 0 

K 



dimensionless 


Subscripts : 
a air 


c compressible 
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i inc ompr e s s ib le 

m maximum 

s surface 

w water 

x horizontal component 

y vertical component 
0 free stream 


ANALYSIS 


The differential equations that describe the motion of a droplet 
with respect to a cylinder have been derived in reference 3 and are pre- 
sented here in the following form: 


dv x 

dT 


C P Re 1 . , 

~ K K “ v x) 


dv^ 

dT 


CpRe i 
24 K 


(Uy - 


V 


(i) 


where 

v = 2 Pw a2u 
9 H 0 L 


( 2 ) 


and the Reynolds number Re can be obtained conveniently in terms of 
the free -stream Reynolds number 


Re 0 



(3) 


such that 



The differential equations (l) state that the motion of a droplet 
is governed by the drag forces imposed on the droplet by virtue of the 
relative velocities between the droplet and the local air in motion along 
the streamlines around the cylinder. The compressible-air flow around a 
right circular cylinder of infinite extent is described by the following 
equations : 
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The first bracketed term describes the incompressible flow around the 

cylinder. The term with the factor and the term with the fac- 

tor are, respectively, the first- and second-order corrections for the 
effect of compressibility on the incompressible flow field. The expres- 
sion for the velocity potential, from which equations (5a) and (5b) were 
derived, was obtained from reference 7. In equation (5), a is measured 
counterclockwise from the x-axis of the cylinder, as is shown in figure 1. 

Equations (l) to (5) are written in dimensionless form in order to 
maintain the number of calculations to a minimum and to simplify the 
presentation of the results. The equations apply to the motion of drop- 
lets in a plane perpendicular to the axis of the cylinder, which is 
located at the origin of the coordinate system (fig. l). At an infinite 
distance ahead of the cylinder the uniform air flow carrying the cloud 
droplets is assumed to be approaching the cylinder from the negative 
x-direction and parallel to the x-axis. All the distances appearing in 
the equations are ratios to the cylinder radius L, “which is assumed to 
be the unit of distance. The velocities appear as fractional parts of 
the free -stream velocity U. Time is expressed in terms of the cylinder 
radius and free-stream velocity such that the unit of time is the time 
required for a droplet to travel a distance L at velocity U. 

The more important assumptions that have been necessary in order to 
solve the problems are: 

(1) At a large distance ahead of the airfoil (free-stream conditions) 
the droplets move with the same velocity as the air. 

(2) The droplets are always spherical and do not change in size. 

(3) No gravitational force acts on the droplets. 

(4) The ratio of the local viscosity p to the free-stream viscosity 
Pq, appearing in equation (4), is unity. This assumption is valid, 

because viscosity is not appreciably affected by the order of magnitude 
of the changes in temperature found along a trajectory. 


METHOD OF CALCULATION 

The nonlinear differential equations of motion (equations (l)) are 
difficult to solve by ordinary means, because the velocity components 
of the air and the term containing the coefficient of drag are variable 
at all points in the flow field and cannot be related in closed form. 
Simultaneous solutions for the two equations were obtained with a mechan- 
ical analog based on the principle of a differential analyzer (refer- 
ence 3). The answers were obtained in the form of plots of the droplet 
trajectories with respect to the cylinder, as shown in figure 2. 
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The second- quadrant section of the cylinder is outlined in fig- 
ure 2. As is shown in the figure, the ordinate scale was expanded 
approximately four times the abscissa scale in the mechanical analog in 
order to gain accuracy in determining the points of impingement of the 
droplets on the cylinder surface. The droplets were assumed to start at 
an infinite distance ahead of the cylinder at an initial ordinate of 
Yq (fig. l) , but only the last 5 radii of distance in the x-direction 

were plotted by the machine. The starting conditions at x = -5 were 
obtained by use of equations described in references 2 and 3. The ref- 
erence angle 9 in figure 2 is measured clockwise from the negative 
x-axis of the cylinder. 


Trajectories, such as those shown in figure 2, computed for several 
combinations of values of K and Rsq permit the evaluation of the 

amount and extent of water-droplet impingement on cylinders. The cylinder 
data are presented herein in terms of the dimensionless parameters K 
and cp used in reference 3. The parameter K (equation (2)) has been 
termed the inertia parameter because its magnitude directly reflects the 
external force required on a. droplet to cause a deviation from the 
original line of motion of the droplet. The dimensionless parameter cp 
is defined 


_ Be Q 2 _ 18LU ( P s^ \ 
9 ~ K P v \ H / 0 


( 6 ) 


The magnitude of 9 is a measure of the deviation from Stokes ' law for 
the forces acting on the water . droplets. Stokes' law. was derived for 
slow translatory motion of a small sphere in an incompressible viscous 
fluid and applies precisely in the limiting value of cp = 0. 

The parameters K and cp together are sufficient to define the 
conditions of a particular impinging droplet. Because of the nature of 
the dimensionless parameters, an infinite number of combinations of 
values of the variables,, such as free-stream velocity, cylinder size, 
droplet size, and others (equations (2), (3), and (6)), would apply to 
the particular value of the dimensionless parameter under consideration, 
A system of equations for the evaluation of dimensionless parameters in 
terms of variables with units commonly employed in aeronautics is pre- 
sented in appendix C of reference 3. 

The compressible -air velocity components in the flow field (equa- 
tions (5)) were calculated at a Mach number of 0.4, because the effect 
of compressibility for completely subsonic flow is considered greatest 
at or near the flight critical Mach number. The value obtained for the 
flight critical Mach number depends on the number of correction terms 
applied to the incompressible flow field and appears to approach 0.4 as 
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a limit. This value of Mach number corresponds to a speed of approxi- 
mately 290 miles per hour at an altitude of 10,000 feet. Results pre- 
sented herein involve the maximum effect of compressibility, which 
restricts the free -stream velocity to a value compatible with a Mach 
number of 0.4. 

Because of the difficulties in making precise corrections to the 
coeffieient-of-drag term in equations (A) and to the density terms in 
equation (4), the problem was solved in two steps. The first step was 
to survey a range of values of the parameters K and cp to determine 
the conditions. most affected by compressibility. For this survey the 
term Cj)Re/24, containing the coefficient of drag for the droplets, 

required in equations (l), was assumed to be a function of the local 
Reynolds number Re only and not related to the Mach number. The values 
for Cj)Re/24 were obtained from tables in reference 2. Also, the den- 
sity ratio appearing in equation (4) was assumed to be -unity for this 
survey. The second step was to apply first-order corrections to the 
coeff icient-of-drag term and to equation (4) for density variations 
along the particular trajectory that showed the largest effect caused 
by the compressible -air flow field. The difference in the trajectories 
caused by the correction was less than the probable error of the analog. 


RESULTS AND DISCUSSION 
Presentation of Results 

The solid lines in figure 2 are the trajectories for droplets in a 
compressible -air flow field calculated at the flight critical Mach number 
for the cylinder (Mq = 0.4), and the dashed lines are the trajectories 

for droplets in an incompressible flow field. The set of lines labeled 
A are trajectories of droplets with large momentum, that is, large in 
size and moving with high speed. The momentum of these droplets is so 
large that the trajectories are nearly straight lines. One combination 
of operating conditions for which the trajectories labeled A would be 
applicable involves a cylinder l/2 inch in diameter moving with a free- 
stream velocity of 290 miles per hour through a cloud composed of uniform 
droplets approximately 26 microns in diameter at a 10,000-foot altitude. 
The trajectories labeled A are tangent to the cylinder at the point of 
interception. All droplets with the same values of parameters K and cp 
having trajectories below these trajectories strike the cylinder; 
whereas, all droplets with the same values of cp and K having trajec- 
tories above these lines will miss the cylinder. The initial ordinate 
y^ of the droplets at x = - oo (fig. l) is given in figure 2. 

A comparison of the trajectories computed for the compressible -air 
flow field with those for an incompressible flow field shows that for 
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droplets with large values of the inertia parameter K (trajectories A, 
fig. 2) the effect of compressibility is imperceptible. The effect of 
compressibility is also negligible on trajectories of droplets with low 
values of K, as shown by the set of tangent trajectories labeled B. 

These trajectories would apply, for example, to a cylinder 1 inch in 
diameter moving with a free -stream velocity of 290 miles per hour through 
a cloud composed of uniform droplets approximately 4 microns in diameter 
at a 10,000-foot altitude. 

The greatest effect of compressibility on tangent trajectories is 
found in the middle range of the droplet inertia parameter K, as shown 
by the set of trajectories labeled C in figure 2. Both trajectories 
are tangent to the cylinder, but the droplets had to start at different 
ordinates in order to arrive tangent. For the droplets having trajec- 
tories labeled C, the collection efficiency (E m = yQ^ m ), defined as the 

ratio of the actual water impinging on the cylinder to the total water 
in the volume swept out of its path by the cylinder, is 0.501 if com- 
pressibility is accounted for, as compared with 0.514 when compressibility 
is ignored. 

The collection efficiency E m for a cylinder in a compressible-air 
flow field at approximately the flight critical Mach number (Mq = 0.4) 

is presented in figure 3 as a function of inertia parameter K and 
dimensionless parameter cp. 


The effect of compressibility on the collection efficiency is shown 
in figure 4 as a percent difference. The difference is obtained from 
the relation 


E 


percent difference = 


■“ii. 
E. 


E. 


'm,c 


m,c 


where the values for E m ^ j_ are obtained from the collection efficiency 
presented in reference 3 for the incompressible flow field.,. The greatest 
difference is about 3 percent at the conditions of K = 5 and 
cp = 50,000. The following tabulation summarizes the effect of compres- 
sibility as a percent difference between results computed in a compres- 
sible and an incompressible flow field. Typical values of dimensions 
and flight conditions are involved in the examples given. All calcula- 
tions were made for 290 miles per hour, which is near the critical speed 
for cylinders at the 10,000-foot NACA standard altitude. 
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Inertia 

parameter 

K 

Cylinder 

diameter 

D 

(in.) 

Droplet 

diameter 

d 

(microns) 

Dimension- 

less 

parameter 

cp 

E . - E_ „ 

tp \ 100 J 

*m,c 

(percent) 

0.5 

1/2 

3 

740 

0.3 

.5 

1 

4 

1,480 

.4 

.5 

3 

7 

4,440 

.5 

.5 

6 

9 

8,880 

.5 

1 

1/2 

4 

740 

.6 

1 

1 

5 

1,480 

.8 

1 

3 

9 

4,440 

.9 

1 

6 

13 

8,880 

1.0 

5 

1/2 

9 

740 

1.5 

5 

1 

12 

- 1,480 

1.9 

5 

3 

21 

4,440 

2.1 

5 

6 

30 

8,880 

2.3 

5 

18 

50 

27,650 

2.8 

5 

30 

70 

44,400 

2.9 

20 

1/2 

17 

740 

.5 

20 

1 

24 

1,480 

.6 

20 

5 

40 

4,440 

.7 

20 

6 

60 

8 ,880 

.8 


Discussion of Results 

From the results of the tabulation and of figure 4, it may be con- 
cluded that, for a cylinder the effect of compressibility of the air on 
the droplet trajectories is not very great up to the flight critical 
Mach number; and for most practical problems involving measurements of 
droplet impingement on cylinders, the compressibility effect may be 
ignored. 

The reason that the effect of compressibility on droplet trajec- 
tories is nearly negligible is found in the comparison of the 
compressible-air flow field with the incompressible flow field. The 
effect of compressibility is most pronounced on those streamlines close 
to the cylinder. The compressibility effect decreases very rapidly with 
distance away from the surface of the cylinder and everywhere becomes 
negligible 1 cylinder radius away from the surface. Furthermore, the 
effect of compressibility is of consequence only over the top portion of 
the cylinder, as shown in figure 5. The comparison between the compres- 
sible and incompressible tangential air -flow velocities on the cylinder 
surface (fig. 5) shows that the effect is pronounced only for 9 >60°. 
The angle is measured from the same reference line as in figure 1. 
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Droplets forming the solid-line trajectory labeled A in figure 2 
are acted on, during the last portion of the trajectory, by air stream- 
lines that have a different pattern from the air streamlines for an 
incompressible flow field j however, the inertia of these droplets is so 
large that the last portion of their path is not influenced appreciably 
by the compressibility of the air flow. The droplets with low inertia 
(trajectories labeled B, fig. 2) are always in the flow-field region 
where there is no appreciable difference between the compressible and 
incompressible air flows. The droplets in the middle range of inertia 
have the last portion of the trajectories in a region where the effect 
of compressibility on the air flow is of some consequence. The inertia 
of these droplets is low enough to be influenced. Droplets with trajec- 
tories labeled D and E have the same middle range of inertia as 
those with trajectories C. Trajectories D and E are not influenced 
by compressibility for the same reason that trajectory B was not 
influenced. 

The sample trajectories presented in figure 2 summarize the effect 
of compressibility over the complete range of values studied for the 
parameters K and <P. The small effect of compressibility is most 
pronounced on collection efficiency (figs. 5 and 4). The effect on the 
maximm angle of impingement 0 m was within the accuracy of the analog. 

The effect on trajectories intermediate between the tangent trajectory 
and the x-axis is imperceptible (trajectories D and E, fig. 2). For 
this reason, the values of the local rate of impingement (the term (3 
in reference 3) are not affected by compressibility. 


Extension of Cylinder Results to Airfoils 

The extension of the results obtained with a cylinder to an airfoil 
is justifiable, because the incompressible flow fields of cylinders and 
airfoils are similarly altered by compressibility. A comparison between 
the compressible and incompressible velocities on the surfaces of two 
symmetrical Joukowski airfoils is shown in figure 6. The curves of fig- 
ure 6 were obtained from results presented in reference 8. The ratio 
of the surface velocity to the free -stream velocity is plotted as a 
function of the chordwise station, given as a ratio to the chord length. 
The results of figure 6 indicate that the effect of compressibility 
becomes less pronounced as the thickness ratio t of the airfoil is 
decreased, and the critical Mach number is increased. 

As was found for the cylinder, the greatest effect of compressibility 
is in the vicinity of the thickest part of the airfoil. Analysis based 
on the work presented in references 5 and 6 indicates that under normal 
icing conditions the impingement occurs on the nose section of the air- 
foil. As with the cylinder, the nose section is not affected to a sig- 
nificant degree by compressibility. Those droplets that impinge on the 
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thick part of the airfoil have such large inertia that small differences 
in air-flow velocity components close to the airfoil do not appreciably 
affect their trajectories. 

Aerodynamic considerations show that the superposition of circulation 
on the air-flow field, as would be caused by placing the airfoil at an 
angle of attack, would tend to cause a shift in the location of the max- 
imum difference between the compressible and incompressible air flows. 

The circulation would shift the maximum difference toward the nose on 
the top surface and toward the tail on the bottom surface, because cir- 
culation increases the air velocity around the top surface, especially 
near the nose, and decreases the velocity near the bottom surface. 

However, a change in angle of attack also causes a shift in the area of 
impingement in the same direction, that is, toward the nose on the top 
surface and toward the tail on the bottom surface; therefore, the area of 
impingement would probably remain in the region where the effect of com- 
pressibility is of no practical consequence. 

The curves of figure 6 are for Joukowski airfoils, which are more 
blunt-nosed than modern high-speed airfoils. The chordwise location of 
the maximum velocity occurs at approximately x = 0.1 for the Joukowski 
airfoil j whereas, the chordwise location of the maximum velocity occurs 
from* x = 0.3 to x = 0.5 on high-speed airfoils, depending on the 
design lift coefficient. The area of droplet impingement, however, is 
at the nose section for a high-speed, airfoil (reference 6), as in the 
case of the Joukowski airfoil (reference 5). The arguments used in con- 
cluding that the effect of the compressibility of air is negligible on 
droplet trajectories impinging on a Joukowski airfoil, for which the 
flow field is a direct transformation from a cylinder, are applicable to 
modern high-speed airfoils, because the air -flow regions with the great- 
est compressibility of the air occur on the airfoil well back of the 
regions of droplet impingement. 


CONCLUSIONS 

The following conclusions are based on the calculation of droplet 
trajectories around a cylinder in both an incompressible flow field and 
a compressible-air flow field at the flight critical Mach number and on 
a study of the air velocities around airfoils: 

1. For a cylinder the effect of compressibility of the air on the 
droplet trajectories is not very great up to the flight critical Mach 
number, and for most practical problems involving measurements of 
droplet impingement on cylinders the compressibility effect may be 
ignored. 
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2. The fact that the effect of compressibility on the impinging 
droplet trajectories is negligible also applies to airfoils. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 11, 1952 
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Figure 2. - Trajectories of droplets impinging on a cylinder. 
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Chordwise station, x, ratio to chord length 


Figure 6. - Comparison of surface velocities on symmetrical Joukovski 
airfoils. Angle of attack, 0°, 
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